Size-independence of statistics for boundary collisions of random walks 
and its implications for spin-polarized gases 
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A bounded random walk exhibits strong correlations between collisions with a boundary. For 
an one-dimensional walk, we obtain the full statistical distribution of the number of such collisions 
in a time t. In the large t limit, the fluctuations in the number of collisions are found to be size- 
independent (independent of the distance between boundaries). This occurs for any inter-boundary 
distance, including less and greater than the mean-free-path, and means that this boundary effect 
does not decay with increasing system-size. As an application, we consider spin-polarized gases, 
such as 3-Helium, in the three-dimensional diffusive regime. The above results mean that the 
depolarizing effect of rare magnetic-impurities in the container walls is orders of magnitude larger 
than a Smoluchowski assumption (to neglect correlations) would imply. This could explain why 
depolarization is so sensitive to the container's treatment with magnetic fields prior to its use. 

PACS numbers: 05.40. Fb 76.60.-k 67.30.ep 



Introduction. Random-walks between (or near) 
boundaries crop-up throughout the mathematical sci- 
ences, from diffusion of particles in a box, to biological 
systems PM], or the gambler's ruin problem [5, 6J. The 
last of these is a first passage problem (the probability 
that the walk hits the boundary at its nth step). For 
a standard one-dimensional (ID) random walk, this first 
passage problem has long since been solved [5J EJ, and 
now most work is for higher dimensions in various ge- 
ometries [7H5], or anomalous walks [1TJHT5]. Other works 
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Figure 1: Sketches of bounded random-walks: (a) a diffusive 
walk and (b) a quasi-ballistic walk. In (a) the boundary- 
collisions are clustered. In (b) the boundary collisions are 
anti-clustered (predominantly equally spaced), because ran- 
dom changes of direction (highlighted with stars) are rare. 



study dynamics of walks within finite or bounded regions 
|15H17j . Here we address a closely related problem; we 
take a standard ID random walk trapped between two 
boundaries (labelled "L" and "R"), and study the statis- 
tical distribution of the number of boundary-collisions in 
a time t. Despite a well known formal connection to the 
recurrence time (first-passage time for a walk starting at 
the boundary) I14| , this problem does not appear to 
have been solved before now. We obtain the full distri- 
bution for arbitrary t, and study the large t limit . 

Our central results, Eqs. (3|4 1, show that the statis- 
tical fluctuations in the number of L and R boundary- 
collisions, TVl and 2Vr, exhibit size-independence for large 
time t. Remarkably, this is a boundary effect which does 
not decay with increasing system size, X, even though 
the average number of such collisions decays like 1/X. 
The X-independences applies for all X, both greater and 
less than the walk's mean-free-path, ^ m f p ; see Fig. [ljt, and 
b. It is a consequence of the correlations between subse- 
quent boundary collisions, with collisions clustering for 
large X and anti-clustering for small X. The cluster- 
ing for large X means the variance of A^l,r is very much 
greater than the average. The only requirement is that 
t is much larger than the mean-free-time r m f p = l mp f/v 
(so the motion is random) and the time to traverse the 
system, i trav (so the walk explores the whole system). 

Below we discuss spin-polarized gases ( 3 He, Xe, etc), 
which are used in a variety of scientific and medical sit- 
uations, and argue that the above size-independent fluc- 
tuations are a crucial source of depolarization. 

Model and Results. A telegrapher equation is a 
standard model for one-dimensional random walkers with 
a finite mean free path, l m f p [5J [T51 HH]- It models a 
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particle moving with velocity v in one-dimension, which 
changes the direction of its motion at random on a 
timescale T m f p = l m f p /v. We assume the walker is be- 
tween two boundaries (see Fig.[lJ, the left (L) boundary 
at x — and the right (R) boundary at x = X. The 
walk reflects whenever it hits the boundary i.e. its ve- 
locity is reversed. Working with dimensionless variables. 
t = t/r m f p and z = x/X, we define P±(z, r) as the prob- 
ability densities that the random-walker is at z and in 
state + or — (i.e. has velocity +v or —v) at time r. The 
telegrapher equation is equivalent to [51 118j. 

&P±(z,t) = Te£P±(z,r) - P ± (z,r) + P T (z,r) , (1) 

where e = VT m f p /X. Eq. ([IJ is a pair of master equations 
for the left or right motion of the walker. The first term 
on the right-hand-side is motion in the direction of travel 
(with velocity e in our dimensionless units), while the last 
two terms account for the changes of direction which oc- 
cur at random (at an average rate equal to one in our 
dimensionless units). The diffusive and quasi-ballistic 
limits are £«1 and e ^> 1, respectively. 

For a walk starting at zq at r = 0, we have P±(z, 0) = 
1(1^(1 — 2a)]S(z—Zo), where a is the probability that the 
initial velocity is positive. We write the boundary condi- 
tions at z — 0, 1 as P+(0, r) = /jl-P-(0, t) and P-(l, t) = 
li r P+(1, t), where /^l and /iR give escape probabilities at 
each boundary collision. In reality /il = /iR = 1, how- 
ever we leave them free so we can use them as counting 
variables to track the boundary collisions |2"CH |2T . The 
survival probability [5] (probability to be in the system 

at time t), is *(t|z ) = j dz[P-(z,r\zo) + P+{z, t\z )] ■ 
The Supplementary Material details the calculation [5] 
of the Laplace transform of the boundary-collision statis- 
tics, with f(s) = / °° dt e~ BT f(r) for any /(r). To get 
the long time behavior, we do as follows. 

1. We show that ^(sjzg) is a generating function for 
the desired quantities (generated by taking deriva- 
tives of tf> with respect to /il,r), and that \&(s|.zo) 
is given by P±(z, s\zq) at z = 0, 1. 



2. We get P±(z,s\z ) at z = 0, 1 by writing Eq. (hj)'s 
Laplace transform as a matrix equation, using a 
Fourier transform in z and a matrix diagonalisa- 
tion. We thereby get the exact algebraic expres- 
sions for (JV l ,r(s|^o))) (Nl :R (s\z )), etc, given in 
Eq. (22-27) of the Supplementary Mxaterial. 

3. For small s, we invert the Laplace transforms, to 
get the statistics for t 3> max[r m f p , itrav]- 

We also find that i t rav — X/v + AT 2 /(u 2 T m f p ). Thus, for 
t 3> max[r m f p , ttrav], the above method gives us 



(AL) = (N R ) = vt/(2X), 
var[JV L ] = var[JV R ] = t/(3r mip ), 
covar[7V L ,iV R ] = -i/(6r mfp ) 



(2) 
(3) 
(4) 



Since the distribution of A'rr is gaussian [5 , all higher 
moments are given by the above second-moments, and 
so exhibit the same AT-independence. Eqs. ([2jj4]) apply 
for any X, from quasi-ballistic {X <C l m f p ) to diffusive 
(X Z m f p ). Although for large X/l m f p the time at which 
the dynamics enters the above long time limit goes like 
^trav ot X 2 . In this diffusive limit, the variance is vastly 
larger than the average, although typical fluctuations are 
of order y / var, and remain much less than the average. 

The index of dispersion is the variance-to-mean ratio 
(VMR), and it tells us about clustering [55]. For the L 
or R boundary collisions, the index is 

varLV LR l 2X f > 1 ( clustered ) forI »|rfP 



31 



1 (Poisson-like) for X — |Z m f p , (5) 
mfp I < 1 (anti-clust.) for X < l mfp 



Fig. [IjL,b show the clustering for X ^> l m jp and anti- 
clustering for X <C / m fp- Remarkably, Eq. (I3J) is exactly 
the same in two such physically different limits. 

We also get the statistical distribution of boundary- 
collisions for arbitrary times, as an algebraic expression 
for its Laplace transform. The Laplace transform of the 
probability density that a walk initially at zq experi- 
ences N + = Nl + N R boundary collisions in a time r, 
is g(N+,s\z ) = (A^ + !)- 1 (d/d^) Ar +$(s|zo)| Ai=0 , where we 

take /il = Mr = M- Since 'I' is a fairly simple function of 
\x, one can evaluate the derivatives for any N + . Then 



g(N+>0,s\z o ) 



E[m(a) +u u m(l-a)](l-u u ) 
~ iVx+i ' W 



o\ N ++ 1 
2\eu v r s 



and g(N + — 0, s\z ) = s (1 — 2m(a) / D ) , where 



u± 



m(a) 



sinh [A] ± Ae 
(s + 1) sinh [A] — Ae cosh [A] ' 
(l + (l-a)s)sinh [A(l - z )] 
+ (1 + as) sinh [Xzq] + aXe cosh [A^o] 
+ (1 — a)Xe cosh [A(l — z )} . 



(7) 



(8) 



for A = ^(2 + s)/e and D = (1+s) sinh[A]+Aecosh[A]. 
Performing the inverse Laplace transform numerically, 
Fig. [2] shows how this distribution evolves in time. For 
t max[T m f p , £trav] h becomes the gaussian distribution 

g(N + ,t\z ) = (27ry t )- 1/2 exp [- (N+ - (N + (t))) 2 /(2V t )] , 
where (N+(t)) and V t = v&r[N + (t)] are in Eqs. 

Intuitive picture of the results. We use the clus- 
tering to explain intuitively the surprising result that 
the variances do not decay at large X/l m f p . We cut a 
long random-walk into many segments each beginning 
and ending at x = X/2 (see Figjlji), each taking a time 
t scg ~ 3ttrav The walk takes a time ttrav ~ X 2 / (vl m { p ) 
to diffuse to a boundary, upon which it recoils to a dis- 
tance Z m f p from that boundary. Then the probability 
that it does not hit the boundary again before returning 
to x = X/2 is about 2l m { p /X. This probability is tiny, 
so the segment contains N c \ ustQr ~ X/(2l m f p ) boundary 
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Figure 2: The evolution of the distribution of boundary col- 
lisions with time, for three different e, given by numerically 
inverting the Laplace transformed distribution in Eqs. ( 6|8 1. 
We plot the distribution of N+, for a particle initially at the 
mid-point between the boundaries, zn = a = 1/2. To make 
the time-evolution of the shape of the distribution clearly vis- 
ible, we multiply the vertical scale by the width of the dis- 
tribution, V"/^ 2 , at each time t, where Vt = vax[N+(t)]. The 
gaussian is that given below Eq. Q for i/ttrav = 10. 



collisions. Thus the statistics are similar to tossing a coin 
every time-period t seg and saying that a "head" is A^ c i uster 
collisions at the L boundary, and a "tail" is A^ c i ustcr col- 
lisions at the R boundary. Then (Nz,) and (ATr) go like 
-^cluster x i/^trav ~ vt / X . However, the variances have 
A'duster in P lace of ^duster, and so go like t/r mip . One 
also has cov&r[N^, ATr] ~ — varfA^]. This simple argu- 
ment gives Eqs. (|2j[4j , except the O[l]-prefactors. 

Difference from Smoluchowski. Eq. ^ is very dif- 
ferent from Smoluchowski 's model of Brownian motion 
[25] . When estimating the number of collisions each liq- 
uid particle makes with a macroscopic object (the particle 
undergoing Brownian motion), he neglected correlations, 
taking var[A^L — A^r] ~ var[A r x J + A^r] ~ (N L + N R ). Yet, 
if each liquid particle performs a random walk, Eq. ^ 
shows that is very far from the truth; since the container 
size, X, is many orders of magnitude larger than Z m f p . In 
fact, Smoluchowski's assumption only works for Brown- 
ian motion due to many-body effects (see below). 

Spin-polarized gases. Such gases, particularly 3 He, 
are used as a spin-filter for neutrons (24] , a precision mag- 
netometer |251 126] or to fundamental spin-dependent in- 
teractions |27l I28j . They are used for magnetic resonance 
imaging in medicine [25] and engineering [30 . 3 He gas 
is typically stored at room temperature and at pressures 
0.1-lbar, in a glass container centimeters across [3T] ; un- 
fortunately, it slowly depolarizes during storage. There is 
a great variety in the quality of the containers; the gas re- 
mains polarized for a few hundred hours in the best con- 
tainers, while it depolarizes in only a few hours in other 
superficially identical ones. The depolarization-processes 
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Figure 3: The random walk of a 3 He atom due to collisions 
with other 3 He atoms. On the left, collisions with the con- 
tainer walls are marked by ellipses. On the right, we show the 
effect of one such boundary-scattering on the atom's spin. 



due to 3 He- 3 He scattering in the gas |31ff3"5] or inhomo- 
geneous external magnetic field [271 EU [SHU] differ little 
between containers. Thus it is likely that magnetic im- 
purities on the container walls |311 14*3 49j arc the origin 
of the huge differences in depolarization times. 

Here we assume that such magnetic impurities (act on 
a shorter range than mean free path) slightly rotating the 
atom's spin each time the atom collides with the walls, 
see Fig. H| Refs. G2B EI] give a microscopic justification 
of this. For a typical container (see above) the atomic 
motion is diffusive with r m f p ~ 10~ 10 s, X/v ~ 10 _5 s, 
and ttrav ~ Is. The spin-dynamics for t of order the de- 
polarization time (typically tens or hundreds of hours), 
is deep in the long time regime t 3> itrav 3> ^mfp- The 
atom's 3D diffusive motion is very well approximated by 
three uncorrelated ID random- walks in the x, y and in- 
directions. We take the spin to be rotated by a random 
angle 0Li at the ith collision with the boundary L, and 
a random angle 9 Ri at the i collision with boundary R. 
These angles are tiny, since each atom only depolarizes 
after very many boundary collisions. We assume that all 
rotation are about the same axis (relaxing this assump- 
tion does not qualitatively change the results [52]). In 
this case, the spin-polarization at time t is 



S(t) 



S cos[@{t)} ) = Soexp [-fvar[8(t)]] , (9) 



where 0(t) = 6 Li + E;=T d m is the total an g le 

that the spin is rotated in a time t. The average in Eq. {9) 
is over all rotation angles at each collision and all possible 
random walks. Averaging over 0Li and 0Rj, and then over 
the number of boundary collisions, N^(t) and Nn(t), we 
arrive at (0(t)) and var[0(t)]. For (0) = (0 L ) = (0 R ) and 
var[0] = var[0 L ] = var[0 R ], we have (0(f)) = (N + (t)) (0) 
and var[0(t)] = (N+(t)) var[0] + v&r[N+{t)} (0) 2 where 
N+(t) = N h (t) + N R (t). Using Eqs. we find the 

polarization decays exponentially at a rate 

wvar[0] (0) 2 



1 



IX 



6r, 



mfp 



(10) 



for weak enough decay that T\ ^> t tlav [53 . The first 
term is a typical boundary effect cx 1/X for container 
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size, X [33H5T]. The second term looks like a bulk effect, 
but is in fact an X -independent boundary effect, originat- 
ing from the X-independence of Eqs. (3]4|. For typical 



3 He cells VT m f p /X < 10 5 , so T 1 1 is vastly more sen- 
sitive to the average spin-rotation at each collision, (9), 
than the spread of the rotations, y / vaxpj. 

Ref. [37] had a similar result to Eq. ( ]l0| f or quadrupole 
fields, using diffusion equations. Ref. [27 \ got a similar 
result for magnetic-fields in a region within A <C X of the 
wall, using diffusion and Redfield approximations . Their 
result was for A 3> Z m f p (for shorter distances the motion 
is ballistic and does not obey a diffusion equation) . They 
were surprised to find that Monte Carlo simulations of 
random walks showed the same behaviour for A < Z m f p as 
for A 3> l m fp ■ The origin of these paradoxical boundary- 
effect (which did not decay with increasing system-size), 
was not clear for A < l m f p . Our above analysis shows rig- 
orously that such size-independent boundary-effects are 
a generic property of random walks. 

Comparison with experiments. Experiments |48[ 
151] showed a strong reduction of Ti when the container 
had previously been placed in a strong magnetic field. 
"Degaussing" |3H] the container returned T\ to its original 
value. This indicates a low density of magnetic impuri- 
ties on the walls of their container, which the strong field 
aligned and degaussing un-aligned. This strong depen- 
dence of the depolarization on the history of the cell rules 
out bulk effects as the dominant source of depolariza- 
tion. However, it fits with our model of boundary effects. 
The alignment of different impurities causes (9) to grow 
((9) = if they are randomly-oriented), and Eq. (10) 
show its extreme sensitivity to the value of this (9) . In 
contrast, the theory in Ref. [48 assumed no correlations 
between an atom's scatterings from the impurities. So it 
does not explain why Ti should depends so strongly on 
the alignment of the field of different impurities. 

Experiments to-date included a uniform magnetic field 
inducing Larmor spin-precession at frequency u 3> i^ av . 
Our above results are for the motional narrowing |54] 
regime, wi trav -C 1; and so only give a qualitative expla- 
nation of the experiments. We have, however, applied 



our method to wi tr av 1, and find that 1/Xi goes like 
the square-root of pressure. This coincides with the pre- 
diction for short-range forces near the walls [57] , and fits 
reasonably well to experiments in Ref. |51j . 

Applicability of Smoluchowski's assumption. 
Smoluchowski [53] applied his assumption (discussed 
above) to the momentum-transfer between atoms doing 
random walks and a macroscopic object (which under- 
goes Brownian motion as a result). The depolarization 
that we discuss can be considered as a spin-transfer be- 
tween atoms and a macroscopic object (the container 
walls). Despite superficial similarities, these two are very 
different. Inter-atomic scattering in a gas rapidly re- 
distributes the momentum of any given atom to other 
atoms, making momentum-transfer a many-body prob- 
lem. Green-Kirkwood-Kubo mean field theory shows 
that many-body effects typically suppress correlations in 
the momentum-transfer on a timescale of order a few 
Tmfp- Thus Smoluchowski's assumption to neglect corre- 
lations is not unreasonable. 

In contrast, the rate for re-distribution of spin from 
any given atom to the other atoms (due to inter-atomic 
scattering) is tiny; if it were the only depolarization pro- 
cess, then Ti > 800 hours at gas pressures ~ lbar. 
Then boundary collision correlations occur on scales up 
to ttrav ~ Is, and so are unaffected by many-body effects. 
Thus these correlation are crucial for the spin-transfer to 
the walls, and thus for the depolarization rate. 

Conclusions. For a bounded random walk, we gave 
the distribution of the number of boundary collisions in 
an arbitrary time. Surprisingly, the long time limit ex- 
hibits boundary effects which do not decay with increas- 
ing system size. This could explain the 3 He depolariza- 
tion rate's extreme sensitivity to details of the physics 
at boundaries. It would be interesting to consider cases 
where the walk itself has strong correlations, e.g. sub- or 
super-diffusive dynamics. 
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SUPPLEMENTARY MATERIAL 

Derivation of A^l.r statistics. Here we give the cal- 
culation we outlined in the "Model and Results" section 
of the paper, using techniques in Ref. [5]. For a proba- 
bility density F(N-^, N R , t\z ), that the walker hits the L 
and R boundaries TVl and N R times in time r, we have 
*(t\zo) = En„n r =oVl>r R F(Nl,N k ,t\z ). Then 



[dV/dfiw] =1 



w_ 



[d 2 */d/4 
[d 2 */dAi L dAt F 



(N w ), 

(N W (N W -1)), 
(N h N R ), (11) 



for W £ {L,R}. To find W and its derivatives, we 
work in Laplace space, defining the Laplace transform 
of any function /(r) as f(s) = L /(r)e~ ST dr. We 
multiply Eq. (1) by e~ st and integrate over t from to 



oo, noting that J™ dre 8T -^P±(z,t\z ) 
0\z ) + sP±(z, s\z ). Defining the vector 



P ± {Z,T 



P(z,s\z ) 



P + (z,s\z ) 
P-(z,s\zo) 



(12) 



we have 
d 



dz 
where 



P(z,s\z ) = M(s)P(z,s\z )+F5(z-z ), (13) 



1(8) 



1 f-1 - S 1 

-1 1 + s 



HC-i)' < i4) 



The boundary conditions are hJ j 'P(0, s\zq) = and 
bp P(l, s\zq) — 0, where T indicates a transpose of 



1 

"Ml 



Mr 

-1 



(15) 



The Laplace transform of the survival probability is 
SS?(s\zo) = Jp dz(P + (z, s\zq) + P-(z, s\zo)). Subtracting 
the first line of the matrix equation in Eq. ( 13 ) from 
the second line, gives P + + P_ = — s^ 1 (s4^ [P + — P_] — 
8(z — zo))- Placing this in the integral makes it easy to 
evaluate, then using the boundary conditions we get 

$(a\zo) = s- 1 [1 - (1 - Ml)£Pl- - (1 - m)£PR+] , (16) 



where we use the shorthand Pl- = P—(z = 0, s\zq) and 
P R+ = P + (z = l JL a\z ). 

Thus, to find 9(s\zo), we do not need P(z, s\zo) for all 
z, we only need its values at the boundaries z = 0, 1. In 
what follows, we refer to these values using the shorthand 
P L ee P(z = 0,s\z Q ) and P R ee P(z = l,s\z ). To 
get these, we define the Fourier transform of f(z, s) as 
f(k,s) = dze ikz f{z,s). We apply this to Eq. H, 



using J Q L dz e ikz ±P{z, s\z ) = e ik P R - P L - ikP{k, s\z ). 

This gives an equation for P(k, s\zq), which we write in 
the basis where M is diagonal, using 



-i i A 

-A 



c + (l + s — s\) -c_| 
c_(l + .s + eA) -c_ 



where ±A are M's eigenvalues, with eA = W s(2 + s), and 
c± = [(1 + s T eA) 2 - l]- 1 ' 2 . The result is 





Figure SI: On the left we show the short-time behaviour of 
var[iV_(r)] (quoted "a"), (N+(t)) ("b") and var [AT+(r)] ("c") 
where N± — N-r ± Nl. We take a = zo = 1/2, for which the 
mean time to first collision (mtfp) with the boundary, r mt f c = 
(2s)" 1 + (2eY 2 (so r mt f c = 30, 0.75, 0.11 for e = 0.1, 1, 5 
repsectively) . On the right we show the variance ratio for 
various e. Dashed lines are the long-time asymptotes given 
by Eqs. (2-4). On the left we improve these asymptotes by 
adding small r-independent terms 0[l,e -2 ] as fit parameters. 







A + ifc 






-A + ifc 
+ e fco VF 



P(k,s\z ) 

+ VP L - e ife VP R , (17) 



Eq. (17 1 is true for all k, so it must be true for k = iA. 



Then P(fc, s\zq) drops out of the upper elements in this 
vector equation (assuming P does not diverge at fc = 
iA), leaving an equation for Pl and Pr. For k = — iA, 
P(fc, s\zo) drops out of the lower elements, giving us a 
second equation for Pl and Pr. These equations are 



v^Pr 



e A v T P R 



T Pi 

T 



r + L 



F. 



e A2 °v T F. 



(18) 
(19) 



where we define the vectors 
s - eA , -1) and v T ee (0, 1)V = c_(l + s + eX , -1) 



(1,0)V 



R4 



Eqs. (18 191 contain four unknowns; Pl+, Pl-, 
Pr_. However the boundary conditions, Eq. (15|), en- 



able us to write Pr = Ur Pr + and Pr = Ul Pl-, where 

UL=f Ai 1 L V (20) 



Ur = 



1 

Mr 
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Interior of container Glass wall of container 

- Cavities in glass 



Magnetic impurities 



Figure S2: A 3 He atom hitting the wall may enter micro- 
scopic cavities in the wall, getting trapped there for some 
time. This trapping may be geometric (bouncing as shown) or 
chemical (van der Waals bonding to the glass). The trapping 
time is random, and uncorrelated between different collisions 
at the same place on the wall (e.g. longer for trajectory 2 than 
1), The spin's rotation 9 is given by the field at the trapping 
site (due to magnetic impurities), and the trapping time. 




and P R +, Pl- are the two quantities we need for Eq. ( 16 I. 
Substituting this into Eqs. ( 18|19 1, we get equations for 
P R + and Pl-> which in matrix form are 



XT T 

e A viu R -v!u L 

\ T T 

e A viu R -viu L 



Pl- 



c -A 20v T F 



(21) 



Inverting this equation we find that 



-,Az 



Pr+ = 



P,_ 



[vTu L ][vT F ]„ e -A Z Q[ v T UL][v T F] 

e A [vT UL ][v T u R ] 



- a [vIu l ][ V t Ur ] 



-A(l-*o)r v T 



v^u R ][v T F] - e 



A(l-zo) 



[V^U R ][ V T F ] 



c A [vT UL ][vZ;u R ] - e- A [vZ!u L ][vTu R ] 

Here v|u R = l+s-^ R =FeA, v|u L = (I+s^l-ITMleA, 
and v±F = e _1 (l + as =p aeX), so we get 

P+(l,s\z ) 



sinh[Az ] + B£ Xe cosh[Azo] 



P-{0,s\z Q ) = 



2eD 

sinh[A — Az ] + B^ Xe cosh[A — Az ] 
2eD 



where it 



(1 - fM W )(s + 2) T (1 - 2o)(l + 
w — (1 + Mw) T (1 - 2a) (1 - fj, w ) and 
D = [(1 - fi R ){l - /i L ) + s(l + /i L /i R )] sinh[A] + (1 - 
MlA'r) Ae cosh[A] . Substituting these into Eq. ( [l6| , we get 
the Laplace transformed survival probability, ^(sjzo). 



(N 2 ± (s\z )) 

/(1± l)Aecosh[A] A . , 



± 



=F 



Ae (a cosh [A(l — z )] + (1 — a) cosh [Xz ]) 



c 2 nj 



ih[A] 



(1 + as) sinh [A(l — zq)] + (s + 1 — as) sinh [Azq] 



sinh[A] 



(23) 



(A>+( S |zo)A>-(s|z )) 

,~ , . xX s sinh [ Al — Xe cosh [Al 
= <*-H«>)) + 2s^ [X] 

x ^(1 — 2a)s(sinh[Az ] — sinh[A(l — z )]) 

-Ae( cosh[Az ] + cosh[A(l - z )])\ (24) 

To write these in terms of the number of L and R bound- 
ary collisions, we use (N R (s)) = |(JV + (s)) + §(iV_(s)), 

(JV L ( S )> = |(iV + ( S ))-i(JV_( S )) ) while 

(^(s)) = |<^(,s)) + i<^( S )> + !(#+(«)#_(*)>, (25) 
= l(Nl(s)) + \(Nl(s)) l(N + (s)N.(s)), (26) 
(N L (s)N R (s)) = l(Nl(s)) l(N 2 _(s)). (27) 

where for compactness we do not show explicitly that all 
these quantities depend on Zq. 

For s <C min[l,e 2 ], one can perform a small s ex- 
pansion (we must keep terms up to s~ 2 ) and perform 
the inverse Laplace transform, getting the results for 
r 3> max[l, T trav ]. We give these results in Eqs. (2- 
4), neglecting i-independent terms of 0[1] and 0[e~ 2 ], 
equivalent to a small shift of time, t — > t + 0[r m f p ,tt rav ]. 
Fig. |Sl| shows the finite-time behavior of these quantities, 
found by numerically inverting the Laplace transform of 
Eq. (22 24 1. They go fairly rapidly to their long-time 
linear-T behavior, with var[A f _(r|zo)]/var[7V + (r|2;o)] — > 3 
for long-times. 

The mean time to first collision (mtfc) is r mt f c — 

J^drT (dP(0, 0, T\z )/dr) = *(0|2!o)Ul,r=0) and " as a 
simple form; r mt f c = Tb a ii + t^s, where the times in the 
ballistic and diffusive limits are Tbaii = [a(l — Zo) + (1 — 
a)zo]/e and Tdiff = Zq(1 — zn)/(e 2 ). The typical time to 
traverse the system, T trav ~ 2r m tf c (zo = 1/2); which is 
i trav ~ X/v + A" 2 /(u 2 T m f p ) in dimensionful units. 

We conclude by pointing out that we never needed 
the solution of Eq. (13). However for completeness we 

point out that it is P(z, s\z ) = e zM(s) P(z = 0, s\z ) + 



fives exactly. It is convenient to define N±(s\z ) 
Nr(s\zq) ± A^l(s|zo), for which we have 



Laplace transforming Eqs. ([Tl| we evaluate the deriva- e z ° M ^F8(z~z ) , where 9{z) is a Heaviside step func- 
tion, and the vector P(z = 0, s\zq) is given above. 

Intermediate times in diffusive regime. In the 
diffusive regime wr m f p <C X, we can look at the intermedi- 
ate t regime, defined by r m f p <C t <C T trav . We follow the 
above derivation, but now consider e 2 <C s <C 1. Then 
the collision statistics are very sensitive to the walker's 
initial position, xq; if xq is too far from a boundary, there 



(N±(s\z )) = 



(1 — 2a)s (sinh [A(l — z )] qp sinh [Xzo]) 

2s 2 sinh[A] 
Ae (cosh [A(l — z )] ± cosh [Xz ] 



2s 2 sinh[A] 



(22) 
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).0 0.5 1.0 1.5 2.0 2.5 3.0 
pressure [bar] 

Figure S3: Data from Fig. 6 of Ref. [50] (squares are sample 
C#4 and triangles are C#13), with ap 1 ^ 2 (our theory), where 
a is chosen for the best fit. We do not include the expected 
crossover to a 1/p behavior at higher p; instead we do not fit 
to the points where this effect may be significant (two points 
marked in gray for sample C#4). The remaining points are as 
well fitted by our ^/p prediction as by the fitting in Ref. [50] 
(a theory introduced in Ref. [48] which is linear in p at small 
p, with the crossover to 1/p at higher p). More experiments 
are necessary to distinguish between the theories. 



are no collisions. Averaging uniformly over xq, we find 
(TVl) and (ATr) are as in Eq. (2), while 

var[7V L ] = var[7V R ] = vt 3 / 2 /(Xr^), (28) 
covar[7V L ,iV R ] ~ - (N L ) 2 = -{vt/Xf. (29) 

Thus they decay with increasing system-size, with the 
covariance having a much smaller magnitude than the 
variances. Unlike in the long time limit, the distribution 
is not gaussian see Fig. 2, and the typical fluctuations, 
rJ, are much larger than the averages, (JVl,r)- 
As a result the variance gives limited information about 
the nature of the distribution. 

Justification of our model for 3 He. In experiments 
[48,50], boundary-induce depolarization is likely due to 
localized magnetic impurities, see Fig. |S2[ However at 
each collision with a given region of the boundary, the 
atom will be trapped there [49] for a different time (see 
Fig. S2), thereby acquiring a different rotation angle. As 



a result, the rotation angle is not directly related to the 
position at which the atom hits the boundary. Thus it 
seems a reasonable simplification to assume all scatter- 
ing at a given boundary by a single distribution of angles 
averaged over that surface. We make this simplification 
in our model, characterizing the distribution by (9) and 



var[6>] . The central message of our discussion of polarized- 
gases was that the depolarization rate is crucially depen- 
dent on the correlations between boundary scatterings. 
At worst, the above argument under-estimates such cor- 
relation effects, by neglecting correlation due by multiple 
scatterings in the same region of the boundary. 

Depolarization with Larmor precession. The 
spin of 3 He in an external magnetic field precesses at 
a rate uj. This was about 10 5 s _1 for the experiments in 
Refs. [48,50] (the other parameters were T m f p ~ 10 _10 s, 



X/v ~ 10 _5 s, and t 



trav 



Is). Thus we want the der 



3po- 

larization rate, for r m f p <C <C ttrav Here Tj~ is 
the rate of relaxation along the axis of the external field 
(the perpendicular decay rate, T^T 1 , is much faster). 

We argue that w™ 1 acts as a cut-off on the correlations 
of the rotation angle. If the spin rotates clockwise about 
the y-axis at every boundary collision, the angle to the 
z-axis (the external field axis) will systematically grow 
for times <C u , however after a time of order (2w) _1 
the spin will have precessed 180° about the z-axis. At 
this point clockwise rotations about the y-axis will reduce 
the angle between the spin and the z-axis. With this cut- 
off at the variance in the rotation angle is the sum 
of the variances acquired in each time-slice of order 
so var[6(t)] - uit [(N + (uT^)) var[6»]+var[Ar + (cj- 1 )] (9) 2 ]. 
We substitute in Eqs. ( 28|29 1, and get the depolarization 
rate for T mfp < uj^ 1 < t tra v, 



(30) 



Hence in this regime, 1/Ti decays with increasing system- 
size, X. Since T~f p oc gas pressure, p, the second term 
goes like yfp. There is no experimental consensus on the 
p-dependence in such systems [48,50]. While we do not 
know all the experimental conditions, the data in Fig. 6 
of Ref. [50] can be fitted with a y/p, see Fig. |S3l 

For typical experime ntal parameter, the preiactor on 
the second term in Eq. (30 1 is of order 10 5 / 2 ~ 300 times 
larger than that of the first term. Compared with the 
case with no external magnetic field (where the prefactor 
on the second term was 10 5 times larger than that of the 
first term) , the magnitude of the second term is greatly 
reduced. This is a consequence of the fact that corre- 
lations are absent on timescales larger than None 
the less, the remaining correlations are sufficient to make 
1/Ti much more sensitive to average angle, (9), than to 
the spread of angles, yjvai[9]. 



